Evolution from BCS to BKT superfluidity in one-dimensional optical lattices 
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We analyze the finite temperature phase diagram of fermion mixtures in one-dimensional optical 
lattices as a function of interaction strength. At low temperatures, the system evolves from an 
anisotropic three-dimensional Bardeen-Cooper-Schrieffer (BCS) superfluid to an effectively two- 
dimensional Berezinskii-Kosterlitz-Thouless (BKT) superfluid as the interaction strength increases. 
We calculate the critical temperature as a function of interaction strength, and identify the region 
where the dimensional crossover occurs for a specified optical lattice potential. Finally, we show 
that the dominant vortex excitations near the critical temperature evolve from multiplane elliptical 
vortex loops in the three-dimensional regime to planar vortex-antivortex pairs in the two-dimensional 
regime, and we propose a detection scheme for these excitations. 
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Ultracold atoms in optical lattices are ideal systems to 
simulate and study novel and exotic condensed matter 
phases. Remarkable success has been achieved experi- 
mentally with Bose atoms loaded into three-dimensional 
(3D) optical lattices, where superfluid and Mott-insulator 
phases have been observed In addition, experimen- 
tal evidence for superfluid and possibly insulating phases 
were found for fermionic atoms ( 6 Li) in 3D optical lat- 
tices Compared with the purely homogeneous or 
harmonically trapped systems, optical lattices offer ad- 
ditional flexibilities and an unprecedented degree of con- 
trol such that their physical properties can be studied 
as a function of onsite atom-atom interactions, tunneling 
amplitudes between adjacent sites, atom filling fractions 
and lattice dimensionality. 

For instance, in strictly two-dimensional (2D) systems 
the superfluid transition for bosons and fermions is of the 
Berezinskii-Kosterlitz-Thouless (BKT) type @, i|. This 
phase is characterized by the existence of bound vortex- 
antivortex pairs below the critical temperature Tbkt, 
and evidence for it was recently reported in nearly 2D 
Bose gases confined to one-dimensional (ID) optical lat- 
tices [5]. Thus, it is very likely that one of the next re- 
search frontiers for experiments with fermions in optical 
lattices is also the investigation of such a transition. 

For bosons or fermions, it is possible to study not only 
3D and 2D superfluids as two separate limits, but also the 
entire evolution from 3D to 2D by tuning nearly continu- 
ously the tunneling amplitudes [g, [7| . However, fermions 
offer the additional advantage that their interactions can 
also be tuned using Feshbach resonances without hav- 
ing to worry about the collapse of the condensate, as it 
is the case for bosons. Furthermore, the phase diagram 
of fermions in optical lattices also shows superfluid-to- 
insulator transitions [H, 0, [1(3] like bosons do. 

Anticipating experiments, we study in this manuscript 
the finite temperature phase diagram of attractive 
fermion mixtures in ID optical lattices, and discuss the 



dimensional crossover from an anisotropic-3D BCS super- 
fluid to an effectively 2D BKT superfluid as a function 
of interaction strength and tunneling parameters. We 
show that vortex excitations near the critical tempera- 
ture change from elliptical multiplane vortex loops in the 
anisotropic-3D BCS regime to planar vortex-antivortex 
pairs in the 2D BKT regime. Finally, we propose an 
experiment for the detection of vortex excitations. 

To describe fermion mixtures in ID optical lattices, we 
start with the Hamiltonian (fi = kg = 1) 



(1) 



k,a 



k,k',q 



where the operator a t a creates a fermion with pseudo- 
spin a which labels cither the type of atoms for unequal 
mass mixtures or the hyperfine state of atoms for equal 
mass mixtures. The operator 6 fc q = o£ +q / 2 j a -k+ q /2 j. 
creates fermion pairs with center of mass momentum q 
and relative momentum 2k, while g > and Tk are the 
strength and symmetry of the attractive interaction be- 
tween fermions, respectively. Here, £k,er = £k,o- — £t<r rep- 
resents the difference between the kinetic energy 



£k,, 



k 2 ± /(2m a ) + 2t z . a [1 - cos(k z a z )] 



(2) 



and the chemical potential fi a , and a z is the lattice spac- 
ing along the z direction. We allow for the fermions to 
have different masses m a and tunneling amplitudes i 2i<T , 
but we confine our analysis to equal population mixtures. 
The saddle-point action for this Hamiltonian is 

S (A*,A ) = /3|A | 2 /ff+(l/A/)E{' 3 (^,i-i?k,2) 

k 

+ ln[(l + * kl i)/2] + ln[(l + * k>2 )/2]}, (3) 

where f3 = l/T is the inverse temperature, M is the 
number of lattice sites along the z direction, E^ s — 
(£k+ + |Ak| 2 ) 1,/2 + 7s£k,- is the quasiparticle energy 



2 



when 71 = 1 and the negative of the quasihole en- 
ergy when 72 = -1. Here, £ k ,± = (£ k , T ± £ k „i)/2, 
^k,s = tanh(/3i? k . s /2), and A k = A r k is the saddle- 
point order parameter. 

The order parameter equation is obtained from the sta- 
tionary condition OSq/OAq — 0, leading to 



the coefficients for low momentum are c,- 



1/5 



KJ/ v „ , 

= (f/Af)^|r k | 2 A' k!+ /(2 J B k , + ), 
k 



(4) 



where X k) ± = (X kjl ± X^ >2 )/2 and E ki ± = (-E k) i ± 
£k,2)/2. We may eliminate g in favor of the binding en- 
ergy eb < of two fermions in the lattice potential via 



i/.9-(i/M)E k |r k | 2 /(e k ,T 

actions with range i?o ~ k 
and zero otherwise, leading to 



- e ki x — €t) . For s-wave inter- 
we take r k = 1 for k < ko 



e b = 4i z ,+ - (2t% + /e ) exp(l/G) - 2e exp(-l/G), (5) 

where ^4 is the area in the (x,y) plane, eo = k 2 /{2m+), 
= (t 2 ,| + i Zj j)/2, G = m + Ag/(4:7r) is the dimension- 
less interaction strength, and m± = 2m^mi/ (m^ ± rojj. 
Notice that two-body bound states in vacuum only exist 
beyond a critical interaction strength G c = 1/ ln(eo/f 2 ,+) 
for finite i z ,+, while they always exist for arbitrarilly 
small G in the 2D limit where £ Z) + — > 0. 

Eq. (j4|) has to be solved self-consistently with the num- 
ber equation N , a = —dSo/{/3dfi a ), leading to 



k 



(1 - 7-*k,-)/2 - &,+A' k , + /(2B ki+ )] . (6) 



Solutions to Eqs. (j4|) and (J6j) constitute an approximate 
description of the system only when amplitude and phase 
fluctuations of the order parameter are small, which is the 
case only at low temperatures, although quantum fluctu- 
ations play a role. However, fluctuations are extremely 
important close to the critical temperature T c . 

The derivation of the fluctuation action is accom- 
plished by writing the order parameter as = 
\A \S qfi + X(q) with AO) = \\(q)\e l9( -i\ where \X(q)\ is the 
amplitude and 8(q) is the phase of the fluctuations. Near 
T c , \Aq\ vanishes, and the fluctuation action reduces to 
S fl (X*,X) = PE q X*(q)L~ 1 (q)X(q) + {(3b/2) J2 q \X(q)\\ 
where the quadratic term is 



f 



f 



£ (?) = 

„ nit 



-Xq/2+k,T + ^q/2-k.J. 



g 2M ^-f C q /2+ k ,T + Cq/2-k,i - ivi 



2 ,(7) 



and the quartic term is b = S k [^k,+/(4£ k , ) — 
pY kj+ /(8^l + )}. Here, X k , CT = tanh(/3£ k>(7 /2), Y k , CT = 
sech 2 (/?^k,o-/2) are thermal factors, and X k .± = (-Xk,T ± 
X k ,i)/2 and F k ,± = (F k , T ± F k ,|)/2. 

The analytic continuation ivi — > + i<5 where 
5 — > 0, and a long wavelength and low frequency 
expansion leads to L^ 1 (q) = a + Y^ijl^ij^j ~ 
duj. The momentum and frequency independent co- 
efficient is a = 1/g - (1/M) £ k X K+ \T*\ 2 /(2& i+ ); 



Ek {P 2 (ilA.1 X ^ Y ^ + 4 .iu .-Vk .H (16&.+) + 

+ .v k . a..^. + ^) - (4 iT - 

4.. "5k/ - Cka)]/(8^,+)}|rk| 2 , where ^ = 9£k, CT /^ 
and = d 2 £] sit < 7 /(dkidkj); and finally the coefficient 
for low frequency is d = lim ta) _ E k ^k,+ [l/(4£ k + ) + 
i7r<5(2£ k:+ — w)/o;]|r k | 2 . Notice that, is diagonal for 



c,5« and cj_ = c x = c y ^ c 



J c&J. Here, is 



s-wave symmetry with — 
leading to L~ x (q) = a + c±q' 2 _ + 
the Kronecker-delta. 

We consider first the strong attraction regime (G ^> 1) 
corresponding to fi a < and \fi a \ ~ |e&|/2 3> t z ,+i where 
L~ x {q) = m + A/ (4n\eb\)[ive — ^s(q) + 2/ig] to lowest 
order of q and vt, with 



^s(q) = q 2 _/{2m B ,±) + 2t B . z [1 - cos(<7 z a z )] 



(8) 



After the reseating ^f(q) = y/m + A/(4:TT\e b \)X(q), the 
quadratic term of Sfi describes non-interacting bosons 
with dispersion wg(q), mass m B ,± = m^ + mj in the 



[x, y) plane, tunneling amplitude is, 



along the z direction, and chemical potential /is = 
/if + /ij — eb- Since the quartic term of Sfi is small, 
the resulting Bose gas is weakly interacting, leading to a 
dominant contribution to the number equation 



N 



fl,(7 = ^ n B [wb (q) - Ms] 



(9) 



which is the same for f and J, fermions. Here, ris(x) = 
\j[eP x —\) is the Bose distribution and JIb = ^b — Vh < 
includes the Hartree shift Vh- 

For G ^> 1, Eq. ([9]) leads to Bose-Einstein condensa- 
tion of tightly bound fcrmion pairs at q = with T c — 
(m + /m B ^)e 2 D/\n( T c/tB,z), where e 2 D = k\ D j{2m + ) 
is a characteristic energy of fermions in 2D. Here, k 2 D 
is a 2D momentum defined through the 2D density 
n 2 D = ^2i)/(2 7r )- We also define an effective 3D den- 
sity n 3D = n 2D /a z where n 3D = /c| D /(37r 2 ), and k 3D is 
the 3D momentum. Notice that, e 2 n = [2/c3Da z /(37r)]e3D 
where e 3 B — k 2 D /(2m+) is a characteristic energy in 3D. 

For fixed t ZlCr , Eq. Q shows that T c is a decreas- 
ing function of G. This is most easily seen for a di- 
lute system where 2£b,z[1 — cos(q z a z )] rs g 2 /(2ms, 2 ), 
such that tob, z = l/(2t5 iZ af) is the effective mass 
along the z direction. In this case, Eq. ([9]) gives 
T c w 0.218[2TO + /(TO 2 3ji m Bi2 ) 1 / 3 ]e3£ ) , which reduces to 
the 3D continuum result T c = 0.2I8ei? of equal mass 
mixtures (llj where m+ = m, rriB,± = m Bt z = 2m 
and = ep. However, T c — > asymptotically when 
ts.z — > or ms.z —> oo, which occurs when the bind- 
ing energy becomes very large (|e(,| ^> y/t z ,-\t z ,\,). This 
limit is clearly unphysical and shows the breakdown of 
the Gaussian theory in ID optical lattices, since the BKT 



3 



transition of tightly bound fermion pairs is not recovered 
in the 2D limit, as can be seen in Fig. [1] 
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FIG. 1: Phase diagram of temperature T (in units of tiD) 
versus G = mAg / (in) for a mixture of equal masses with s- 
wave interactions, interaction range fco ~ 10 4 A;2o, tunneling 
t z 0.043e2D, lattice spacing a z m 0.43/xm and planar den- 
sity ii2D ~ 2.5 x 10 7 cm~ 2 , such that fc2zja z ~ 2.2. To is the 
saddle-point or pairing temperature scale. 

To recover the BKT physics in the G ^> 1 limit 
where the paired fermions live in 2D planes, we 
return to the derivation of the fluctuation action 
Sfi with t zt7 = 0. Taking the order parameter 



as 



\rj(x)\]e ie ( x \ where \rj{x)\ 



cor- 



responds to the amplitude fluctuations and 0(x) 
is the phase of the order parameter such that 
| An | ~> |»y(a;)|, we obtain the phase-only action Sfi(0) = 
OS/2) E 9 [«o(T)«f + QiPij(T)qj] 0(q)9(-q). Here, 
the coefficient k q (T) = (1/4) £ k [| A k | 2 # k>+ /£ ki+ + 
+3^k,+/(2i?k +)] is the atomic compressibility, where 
A,± = ± ^k, 2 )/2 with 3^ s = sech 2 (/?i^ s /2); 

while the coefficient of the gradient term 



4m i 



Ekikj 
16Al 

k 



4|Ai 



E 3 



/%,+ 

2^, 



/5 



y k ,i y k , 2 



4^ k .- 



Ek,- 



10) 



is the phase stiffness, where Pij(T) = po(T)Sij for the 
s-wave symmetry. 

This leads to the BKT transition temperature [H, 0] 



BKT 



kPo(T B kt)/'Z, 



(11) 



which needs to be solved self-consistently with Eqs. (j4|) 
and (J6j) in order to determine Tbkt, |Arj| and /i<t as a 
function of G. tn the weak attraction regime (G < 1), 
Tbkt increases with G as Tbkt — (e 7 / 7r )(2| e bk2r>) 1/ ' 2 , 
where 7 ~ 0.577 is the Euler's constant and e;, = 
— 2en exp(— 1/G) is the binding energy in 2D. While in 
the strong attraction regime (G 3> 1), Tbkt saturates 
to Tbkt = (e2£>/8)[l - 2m 2 _/(3m 2 _)]. For equal mass 
mixtures (m_ — > 00), this reduces to Tbkt = 0.125e^ 



which can be seen in Fig. [TJ Here, ep = £21? is the 2D 
Fermi energy. 

To estimate when G induces the crossover from 
anisotropic-3D to 2D behavior, we compare the critical 
temperature T c ^Q auss obtained from the Gaussian theory 
with the critical temperature Tbkt for the BKT transi- 
tion in the strict 2D limit. When G>1, the condition 
Tc.Gauss = Tbkt leads to 

3 



tz,c,ltz,c,l 



C 2 (3/2) 



>12, 



1 - 



2m+ 
3^ 



C2D £6 



(12) 

where C( x ) is the zeta function. This relation reduces to 
t z , c = [C(3/2)/(32 v / 7r)] (e 2 Z5|e b |) 1/2 for mixtures of equal 
mass fermions and equal tunneling. 

We can also relate t z ,a to the depth Vb, CT of the ID op- 
tical latttice potential V a (v) — Vo,<r sin 2 (nz/a z ) where 
t z>a = (2E 1%a /^/n)(a a ) 3 ^ 4 exp(-2^/a^). Here, a a = 
Vo,a/Er,(x where E r ^ = n 2 / (2m cr a1) is the recoil energy. 
In Fig. [21 we show the characteristic t ZlC = t ZjCiCr and 
Vo.c = Vb.c,cr lines which separate the anisotropic-3D from 
the 2D regime, for mixtures of equal mass, equal tunnel- 
ing and s-wave interactions. When G is fixed, the 2D 
regime may be reached from the anisotropic-3D regime 



While 



with increasing Vq — Vo,a ° r decreasing t z 
for fixed Vq = Vb lCT or t z — i Z)tr , the 2D regime may be 
reached from the anisotropic-3D regime by increasing G. 
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FIG. 2: Characteristic (a) tunneling amplitude t z , c (in units of 
£2d); and (b) optical lattice depth Vb, c (in units of E r ) versus 
G = mAg /(An) showing the anisotropic-3D to 2D crossover. 
The parameters are the same as in Fig. [TJ 

Further insight into the dimensional crossover is gained 
by rewriting Sfi with t za 7^ in real space and 
time such that 5// (A*, A) = (l/V) J dtdr ± dzC f i(X* , A), 
where C fl (X*,X) = X*(x)(0 - c z d 2 z )X(x) + 6|A(a;)| 4 /2 
is the Lagrangian. Here, V is the volume, X(x) = 
X(r±,z,t) is the fluctuation field and O — a — 
c±\7 2 L — iddt- Upon discretization z = na z , Sfi reduces 
to the Lawrence-Doniach (LD) action <Sld(A* , A„) = 
[1/(MA)] £„ / dtdr ± C LD (X* n , A„), where 

C LD (X* n ,Xn) = KOX n + %|A„ +1 - A„| 2 + ^|A„| 4 (13) 



is the LD Lagrangian [Tz 



Here, the local field 



A„ = X(r±,z = na z ,t) describes the order parame- 
ter in each plane labeled by index n. Writing a — 



4 



o e(T) with e(T) = (T - T c )/T c , scaling the field 
'0n = \fb/a,o\ n , and defining the correlation lengths 
£q j_ = c±/a,Q and £g 2 = c z /a$, and the charac- 
teristic time To = —d/ao leads to the scaled action 
C LD {r n ^n) = r n iTod t ipn + e{T)\yj n \ 2 + Q,A^M 2 + 
£q z \ip n+ i — i() n \ 2 / a1 + \ip n \ 4 /2, which describes the system 

near T c . Here, C LD (ip*,ip n ) = {b/al)C LD (\* n , A„). Fur- 
thermore, taking ip n = \ip n \ exp(i#„) in the LD action, 
such that \tp n \ — <Pa is independent of position and time, 
leads to the phase-only anisotropic-3D XY model with 
the dimensionless Hamiltonian 

H X Y(r±,n) = K±\a ± \7 ± 6 n \ 2 - 2K z cos(6 n+1 ~ 9 n ) + C, 

where K± = {£,o,±4>o/a±) 2 , K z = (Co^0o/a z ) 2 and 
a± ~ , and C is a constant. The dimension- full 
Hamiltonian is HxY(r±,n) = (b/a 2 t )HxY(r±,n). 

This can be mapped into the vortex-loop representa- 
tion [l3| yielding the dual dimensionless Hamiltonian 

H D = 7T [KM*) ■ Mr') + K ± J z (r) ■ J z {r')]U(R), 

where t/(R = r — r') plays the role of an interaction po- 
tential for the vortex loop field J(r) = [Jj_(r), J z (r)], and 
satisfies the differential equation (Vj_ + ?7 _2 <9 2 )[/(R) = 
— 47T(5(R). Here, r\ = \/K~[/K z is the anisotropy ratio 
and <5(a;) is the delta function. The dual transformation 
maps closed supercurrent flows associated with the gra- 
dients of the phase into the vortex- loop vector J(r), in 
the same way that the electric current flowing on a ring 
can be mapped into a magnetic field vector with the help 
of the Biot-Savart law. 

For large magnitudes of R = (R±,R Z ), the vortex- 
loop interaction behaves asymptotically as J7(R) ~ 

1 / y/[R±/ '(^ajj] 2 + (R z /a z ) 2 , and leads to equipotentials 
in the shape of ellipsoids [R±/(rja±_)] 2 + (R z /a z ) 2 = Uq 2 
when U(R) = Uq. Elliptical vortex loops correspond- 
ing to a nearly toroidal arrangement of the supercurrent 
flow are the large scale excitations formed by a contin- 
uous closed line having the same potential between seg- 
ments with r = — r'. When r\ — > oo, the planes along 
the z direction decouple (2D BKT regime) and the vor- 
tex loops reduce to planar vortex-antivortex pairs. For 

2 < r\ < oo, the system is still nearly 2D, and the dom- 
inant excitations are square vortex loops coupling two 
consecutive planes and planar vortex loops. However, 
in the anisotropic-3D regime when rj < 2, the dominant 
excitations become multiplane elliptical vortex loops. 

In the strong attraction regime (G ^> 1), 77 « 
(a z /a±)y // mB,z/'mB,± 3> 1, and the 2D BKT limit is 
recovered since ms,z 3> ™_b._l- For fermion mixtures 
of equal masses and equal tunnelings, we can rewrite 
this condition as 77 fa y / e^D|efc|/(47rf z ) ^> 1. Since the 
anisotropic-3D to 2D crossover occurs for 77 2, this 



condition leads to t ZjC « (l/87r)(e2_D|eb|) 1 / 2 , which is es- 
sentially the same result obtained by equating the Gaus- 
sian and BKT critical temperatures. 

Vortex-antivortex pairs have been detected in ultra- 
cold atoms using absorption images of an expanding 
cloud [14j , and we expect that similar techniques can be 
used to detect vortex loops as the system evolves from 
anisotropic-3D to 2D regime. Vortex loops should appear 
as dark rings in the image since there are no atoms to ab- 
sorb light in their cores. At temperature T, the ratio of 
characteristic in situ core size of vortex loops in the (x, y) 
plane £±(T) = £0 j_ l e C^) |~ 2 ^ 3 and along the z direction 
UT) = to,MT)\- 2/3 is U(T)/UT) * 0.91 for 77 = 2 
and k2DO-z = 2.2. Since typical values of £o,± ~ 0.5/zm., 
then £_l(T) ~ 2.3/im and £, Z (T) sa 2.5/j,m at tempera- 
tures T — 0.9T C , and vortex loops extend to nearly six 
planes for an optical lattice with a z w 0.43/jm. Smaller 
values of rj or larger values of A^u^z enlarge £ Z (T). For 
parameters 77 = 1.7, ^Diz = 5.0 and T = 0.9T C , the 
ratio £±(T)/£ Z (T) « 0.34 and &(T) w 6.87ml, such that 
vortex loops extend to nearly sixteen planes in optical 
lattices with a z rs 0.43/jm. 

We analyzed the finite temperature phase diagram of 
attractive fermion mixtures in ID optical lattices. At 
low temperatures, we found that a dimensional crossover 
from an anisotropic-3D (BCS) to an effectively 2D (BKT) 
superfluid occurs as a function of attraction strength 
eventhough the tunneling amplitude is fixed. In addition, 
we discussed that vortex excitations change from ellipti- 
cal multiplane vortex loops in the anisotropic-3D regime 
to planar vortex-antivortex pairs in the 2D regime, and 
suggested an experiment to detect their presence. 

We thank the NSF (DMR-0709584) for support. 
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